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Abstract. We explore the existence of time reparametrization symmetry in p-spin 
models. Using the Martin-Siggia-Rose generating functional, we analytically probe 
the long-time dynamics. We perform a renormalization group analysis where we 
systematically integrate over short timescale fluctuations. We find three families of 
stable fixed points and study the symmetry of those fixed points with respect to 
time reparametrizations. One of those families is composed entirely of symmetric 
fixed points, which are associated with the low temperature dynamics. The other 
two families are composed entirely of non-symmetric fixed points. One of these two 
non-symmetric families corresponds to the high temperature dynamics. 

Time reparametrization symmetry is a continuous symmetry that is spontaneously 
broken in the glass state and we argue that this gives rise to the presence of Goldstone 
modes. We expect the Goldstone modes to determine the properties of fluctuations in 
the glass state, in particular predicting the presence of dynamical heterogeneity. 



PACS numbers: 64.70. Q-, 61.20.Lc, 61.43.Fs 
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1. Introduction 

Very slow dynamics is an essential feature of glasses p. In both structural glasses 
and spin glasses slow dynamics is manifested through a dramatic increase in relaxation 
times. This slowdown has been captured in mean field theories, such as the mode 
coupling theory for supercooled liquids [2] and the dynamical theory for mean field spin 
glass models [31 HJ |5l |6] . Even though mean field theories are useful in describing some 
aspects of glassy dynamics, they do not completely capture phenomena associated with 
fluctuations. Fluctuations have been shown, particularly with the discovery of dynamical 
heterogeneities, to be central to an understanding of glassy dynamics [7]. 

Dynamical heterogeneities - mesoscopic regions that evolve differently from each 
other as well as from the bulk - have been found in experimental studies of materials 
close to the glass transition [HI |9l [10] and in simulations of both spin glasses and 
structural glasses [11], |12l [13]. Their presence has been directly observed at the 
microscopic level in experiments on colloidal glasses [TO] and granular systems |14j . 
Understanding the onset of heterogeneities without an apparent structural trigger is 
believed to be key to an understanding of the glass transition |7]. There have been 
several theoretical attempts to explain the emergence of heterogeneous dynamics as 
the glass transition is approached. One of them is a geometrical picture, according 
to which dynamical heterogeneities result from non-trivial structure in the space of 
trajectories due to dynamical constraints [15]. Another proposed explanation is provided 
by the Random First Order transition (RFOT) approach, in which a liquid freezes 
into a mosaic of aperiodic crystals [16]. Here we will explore a different theoretical 
avenue to explain dynamical heterogeneities, which is based on time reparametrization 
symmetry [HI [TBI [19] . 

Time reparametrization symmetry (TRS), the invariance under transformations 
of the time variable t — )■ h{t), was discovered some years ago in the mean- field non- 
equilibrium dynamics of the Sherrington-Kirkpatrick model and the p-spin model [5l|6]. 
The symmetry, which was shown to be present in the long-time limit of the mean field 
evolution, implies that the asymptotic equations do not have a unique solution [5l [6l [20] . 
In more recent studies, TRS has been proved to be present in the long time 
dynamics of the glass state in a short range spin glass model, the Edwards-Anderson 
model [T71 [TH [19]. In this last case, the proof of the symmetry is at the level of the 
generating functional, including all fluctuations. Using the renormalization group (RG), 
it was shown that the stable fixed point of the generating functional corresponding to 
glassy dynamics is invariant under reparametrizations of the time variable. However 
not all models of interacting spins under Langevin dynamics show this behavior. For 
example, in a study of the 0(N) model it was shown that the symmetry is not present, 
even for the long time limit of the low temperature dynamics [25]. The explanation for 
dynamical heterogeneities from TRS is derived from the fact that TRS is spontaneously 
broken by the correlations and responses in the glass state. A spontaneously broken 
continuous symmetry is expected to give rise to Goldstone modes, and these modes are 
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associated with spatially correlated fluctuations of the time variable, which give rise 
to heterogeneous dynamics. The proposal that dynamical heterogeneities originate in 
reparametrization fluctuations of the time variable is supported by positive evidence 
from numerical studies in spin glasses [18] and in structural glasses [231 121] • 

In the present work we go beyond mean field theory and use a renormalization 
group procedure to study the presence of time reparametrization symmetry in the 
long time dynamics of p-spin models with arbitrary interaction range, including all 
fluctuations. We consider a system of soft spins on a lattice, with p-spin interactions. 
The spin couplings are assumed to be uncorrelated Gaussian random variables with 
zero mean. We assume Langevin dynamics for the spins with a white noise term that 
represents the coupling of the spins to a heat reservoir. We set up the calculation 
by writing the generating functional of the spin correlations and responses using the 
Martin-Siggia-Rose (MSR) approach and introduce two-time fields that are associated 
with the spin correlations and responses. To study the long time dynamics we start 
the renormalization group procedure by introducing a short time cutoff tq. We 
systematically increase the short time cutoff by integrating over the two-time fields 
associated with the shortest time differences, thus following a procedure analogous to 
Wilson's approach to the RG. In our case, however, we integrate over fluctuations 
that are fast in time^ not in space. We find three families of stable fixed points. 
The first family corresponds to fixed point actions containing the coupling to the 
thermal bath but not the spin-spin interactions. The fixed points in this family are 
not time reparametrization invariant, and we believe that this family is associated with 
the high temperature dynamics. A second family of stable fixed points that are not 
time reparametrization invariant corresponds to fixed point actions containing both the 
coupling to the thermal bath and the spin-spin interactions. For the third family, the 
spin-spin interaction term is marginal but the coupling to the thermal bath is irrelevant. 
The fixed points in this last family are time reparametrization invariant, and we believe 
that they represent the low temperature glassy dynamics of the model. After obtaining 
these results, we discuss their connection with dynamical heterogeneity in the p-spin 
model, and we speculate on how a similar procedure may be applied to models of 
structural glasses, which have been shown to be connected to the p-spin model [211 ES] . 

The rest of the paper is organized as follows: in Sec. |2] we start by giving 
a description of the model and an illustration of how we derive the Martin-Siggia- 
Rose generating functional; in Sec. |3] we show how we use Wilson's approach to the 
renormalization group to get stable fixed points; in Sec. IHwe study the stable fixed point 
generating functionals and determine which ones are invariant under reparametrizations 
of the time variable; and in Sec. |5]we end with a discussion of our results and conclusions. 
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2. Model and MSR generating functional 



The p-spin Hamiltonian is given by 

ii...ip 

where the ... are soft spins subject to the spherical constraint = ^5 

and the couplings are assumed to be uncorrelated, Gaussian distributed, zero mean 
random variables, P{J} = nii<...<ip y=== exp[- .y4ir,j...ij. The dynamics is 

given by the Langevin equation 

where the {r]i{t)}i=i^,„^N are assumed to be zero mean gaussian random variables with 
the correlation {r]i{t)r]j{t')) = 2T6ij6{t — t') that couple the spins to a thermal bath at 
temperature T. Then the Langevin equation can be written as 



P' h...ip^i 



(3) 



From the Langevin equation we use the Martin-Siggia-Rose formalism [26] and 
write down the noise averaged generating functional 

{Z[{h},{h,}]) = J D4>D^D^DN exp U[^,4>] +J2 dtUt)(P,{t) 



+ihi{t)(t)i{t)] + i ^^[Mto) dtN{t) 

i=l 



i=l 
N 



.1=1 



(4) 



where the I and h are sources and the last two terms in the exponent are due to the 
initial condition and spherical constraint, respectively. The action 0] is given by 



1=1 



to 



V 



(5) 



We now average the generating functional over the disorder in the system. The 
action contains only one term with an explicit dependence on the disorder, which we 
call Lj[0, 0], 



II /, dtJi^,„i^(j)i^{t)(j)i^{t)---(j)i^{t). 



to 



(6) 



We now compute the part of the generating functional affected by disorder averaging, 
exp(Lj[0,0]) = / DJP{J}exp[Lj], (7) 
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with DJ = n dJh...ip- Therefore, after integrating over the disorder we have 

ii<---<ip 



exp(Lj[0,0]) given by 



exp(Lj[0,0]) = exp 



dtdt' 



to J to 



ll...lp 

Oi^.a' €{0,1} ii...ipr=l 



(8) 



(9) 



where we have re- labeled the fields using the definitions = and 4>]{t) = 4>i{t). 

p p 

The constrained variables C and C are given by C = (1 — ctj^), C" = (1 — a-^), 

and the constraints C = C" = 1 enforce the condition that for each of the two times t 
and t', there is a product of fields, of which only one is a and all others are (f) fields. 
We are also interested in introducing two-time fields Q"'" (^1,^2), physically associated 
with two-time correlations and responses. In order to do this we write the number 
one in terms of an integral of a product of delta functions that enforce the condition 

Qr'ih,t,)^cpnti)<Pfit2), 

i,ti,t2 ai,a[€{0,l} ^ 

By writing the delta function in exponential form we get 

1 = / DQDQexp ji^ E / / dt^dt^Qf'KhM) 

x(gr'"^(ti,t2)-e(ti)0rk^2))}, (10) 

where we have introduced the auxiliary two-time fields Qi* ^^2) and the notation 
= 1, 1 = 0. We now obtain the noise and disorder averaged generating functional 

{Z[{k},{hi}]) = j DQDQD(tPD(t)^DifDNe^^{S), 

S — Si -\- Sj + Sgpin + Sgxt ~l~ Sbc + Ssc- (11) 
Here we have written the different terms of the action separately: 

5l[g,40^0']=^E E rr 



dtidt2QT''^'{tiM) 



X (gr"(^i,^2)-e(^i)0r(^2)), 

C=1,C'=1 p 



ll...lp 



dtidt^ 



to 'I to 



E n^: 



(12) 

•(^1,^2), (13) 



ai,,a^^e{0,l}r=l 



SsM4>\ 0'] = - ^ E / dtcfrit) \^~'dt4>\{t) + loom] --,Y.K' 

i=l "^*o P- h...ip 



i=l ■'^0 

ftf ft 



'■11. ..ip 



X r r dtidt,9{ti -t2) E n c(^i)</'r>(^2)(i4) 

ai,,a;.„ G{0,l}r=l 
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rtf 



5e.t[0°, 0'; l,h]^ f'dt \hit)<l>iit) + Mt)(t>]{t) 

Jto L 

N 



Ssc[<P\N;N] = i dtNit) 



N 



i=l 



(15) 

(16) 

(17) 



and we have F = Fq, 700 = —iT and g{t — t') =0 at the start of the RG flow. 

3. Renormalization group analysis 
3.1. Renormalization Group Transformation 

We perform a renormahzation group analysis on the time variables. For simplicity we 
take to — and tf — 00 from now on. We focus on the two-time fields. First, we 
introduce a cutoff in the integration of two-time fields, Tq < 1^1—^21- We then write the 
terms of the action affected by the cutoff: 



Si[Q, Q] = ^ E / dt^dt2 E Or'^H^i, ^2) 



* 0<t^,t2<o 
^()<l'l-'2 



X {Qr '{h,t2)-<l>nhW{h)), 

C=1,C'=1 p 



(18) 



P 



ti...ip 



dtidt2 E UQC '^{tiM- (19) 



0<ti,t2<oo 
T"0<l*l-*2l 



"i.-<e{0,l} r=l 



We define fast and slow fields respectively by Q>V"'(ii,^2) = Q"' "'(^i) ^2), for tq < 



\ti — t2\ < bro and Q<V"'(^i, ^2) = ^"""'(^i) ^2), for fero < \ti — t2\, with b > 1. This 
separation of fast and slow parts of the fields results in a separation in the terms: 

s,[Q, g, 0°, = ^i[Q>, g>, 0^ + 5i[g<, g<, 0°, 0^], (20) 

Sj[Q] =Sj[Q>] + Sj[Q^]. (21) 

Next we calculate the integral /> over fast fields. To do this we use the fact that there 
are no cross-terms between fast and slow fields in the integral: 



J DQ^DQ 



> exp < 



^E 



dtidt2 E Q>i"'{ti,'t2) 



* •ro<|ti-t2l<6To 
0<tl,t2<°° 



x(g>rk^i,^2)-0r(^i)0rki2) 



p2 /. 

E -^n...ip / 
h...ip ... ■'. 



C=l,C'=l p 

dt,dt2 E n Q>i 

0<tl,t2<°° 



■(tl,t2 



(22) 
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Calculating the integral over the Q> fields constitutes undoing the delta function integral 
transformation we used to introduce the two-time fields for the fast modes. Hence, 



/> = exp <j I dtidt2 ^ E n €"(^1)0^ (^2) \ .(23) 

I P- j J / ,1 „_1 



C=1,C' = 1 p 



^0<l*l-*2l<'''r-0 
0<ti ,t2<oo 



/ g{0,l}r=l 



Next we re-scale all the one-time and two-time fields, thus restoring the cutoff to its 
original value Tq 





(24) 




(25) 


W = t, 


(26) 


^rm = b'-^<PTit'), 


(27) 


k{bt') = b^%{t'), 


(28) 


h,{bt') = b^^h'^it'), 


(29) 




(30) 



From the definition of the two time fields in Eq. and the transformations of the fields 
we have 

K„a[ = K,+K'^- (31) 

By rescaling the fields in the part of the action arising from the disorder average {Sj) 
we get 

2 

^j[Q]= - V' E ^n.... / dt[dt', 

X """e"' E f{QT^'''{t\A)- (32) 

"ir.ai^e{o,i} n-v ''=1 

Using the relation between \a,a'i and Aq,' given by Eq. fl3Tl) together with the 
constraints C = C" = 1 we get 

Aj = 2Ao + 2(j9 - l)Ai + 2. (33) 

The other important term we need to consider is Sspin- When we rescale the fields 
we obtain 

N 



-P-Y K 



«l...lp 



'T-%4'l{t') + b'-ioo<P"dt')] (34) 

C=l,C'=l 



dt[dt',b'-'git[ - 1',) y: n </'Tr (^'i)'/>'S'- it'. 

«ir><e{o,i}^'=i 
The result is the following set of fiow equations 

p-i _^ r'-i ^ Y-^b^--', 

700 7oo = loob^^, 

giti - t2) ^ g'{t[ - 4) = 6^'^ [g{b{t[ - t'^)) + Cr,<ibt[-bt',\<b 



(35) 
(36) 
(37) 
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where C-p is defined hy C-p = 1 if V is true and C-p = if P is not true, and 

Kel = Ao + Ai, (38) 
At = 1 + 2Ao. (39) 



If we let b = e'^^ = 1 + dl then the flow equations for F and 700 can be written 
dT 
'dl 

d'joo 



^ = -A„e/r, (40) 
dl 



At7oo- (41) 
dl 

The terms that are left to consider are the constraint terms 5*1 and Ssc^ the 
boundary condition term Sbc, and the coupling to the sources Sext- Physically we 
expect that the constraints represented by 5*1 and Ssc will still be valid in the long 
time limit, and therefore those terms should be marginal. This leads to the equations 
A_^ = — 2 — Xai,a[ = — 2 — Xa^ — X^'^ and Aat = — 1 respectively. It is not obvious 
from physical considerations alone what the behavior of the boundary condition and 
the coupling to the sources should be. The exponents A/, Xh, X^ control their scaling 
behavior, and could in principle be chosen to make the terms marginal, but they do not 
play any role in what follows. 

In order to determine the nature of the fixed points, we need to choose values for 
the scaling exponents Aq and Ai. 



3.2. Choice of Scaling Exponents 

In traditional RG calculations one determines an engineering dimension for the field 
variable by requiring that the free theory action be marginal. In our case none of the 
terms in the action is of the same form as the gradient squared term that is usually 
considered to be the unperturbed part of the action and assumed to be marginal. The 
only systematic way to proceed is to consider what happens when each of the terms in 
the action is marginal. We start by noting that the presence of the spherical constraint 
E^=i[Mt)? = N implies that there is an upper bound on the correlation function 
C{t,t') ~ Q^^{t,t'), and therefore we must have the constraint Ai < 0. The case of 
Ai = corresponds to freezing and the strict inequality corresponds to a decaying 
correlation. The terms in the action that are of interest for our analysis are the three 
terms contained in Sgpin'- the spin-spin interaction, the term containing a time derivative 
and the term coupling the system to the thermal bath. As indicated in Eqs. fl33l) . fl38l) 
and ( 139|) . those three terms have the scaling exponents Aj = 2(1 + Aq + (p — l)Ai), 
Xvei = Aq + Ai and A^ = 1 + 2Ao, respectively. By considering the cases in which only 
one of the terms is marginal we get the results summarized in Fig. [H 

In the case in which the coupling to the thermal bath is marginal, we have a line 
Aq = — 1/2 in the (Ai, Aq) plane. Considering the constraint Ai < and the values of Aj 
and Xjjei we find that there is an interval on this line, Ai < 2{p-i) ' which both the spin- 
spin interactions and the time derivative term are irrelevant. Since the coupling to the 
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Scaling Exponents in the p-spin model for p=3 



Key 



Figure 1. The figure shows the different lines along which each one of the three terms 
in Sapin is marginal for p = 3. The red line corresponds to a marginal coupling to 
the thermal bath (At = 0), the black line corresponds to a marginal time-derivative 
term (A^e; = 0), and the blue line corresponds to a marginal spin-spin interaction term 
(A J = 0). 



thermal bath is marginal then we have a family of stable high temperature fixed points. 
Second, we consider the case where the time derivative term is marginal, corresponding 
to the line Aq = — Ai in the (Ai, Aq) plane. Since we have Ai < 0, the exponent At of 
Eq. (I39l) is always positive, i.e. the coupling to the thermal bath is always a relevant 
perturbation. Thus the fixed points that contain only the time derivative term are 
always unstable. We then consider the case where the spin-spin term is marginal. This 
happens on the line described by Aq = — 1 — (j5 — l)Ai. In the interval 2{^-i) < -^i — 
the time derivative term and coupling to the thermal bath are irrelevant. This gives rise 
to a family of stable low temperature fixed points. Finally, there is the special point 
Ai = 2^}-i) ' -^0 = —1/2, for which both the coupling to the thermal bath and the spin- 
spin interaction are marginal, but the time derivative term is irrelevant, thus allowing 
for an additional family of stable fixed points. 

The above analysis shows that there is a subset of the (Ai, Aq) plane for which 
a high temperature dynamical fixed point family is present. The effective generating 
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functional for this fixed point family is 



Zfp[l,h;T] = {Z[{k},{h^}])^^ = J DQDQD^^D^^DipDN 

xe^phj2[r dt,dt2 E QT'^'Htu t2) h) - <t>T{ti)<t>f{t2) 



-TT.rdt{<pmf 

fOO 

+ dtUt)<f>',it) + iKit)<i>l{t)] 

Jo 

+i E H^I'lito) + i dtN{t) 



N 



.i=l 



(42) 



There is another subset of the (Ai, Aq) plane for which a low temperature interaction- 
dominated fixed point family is present. The effective generating functional for this 
family of fixed points is 

Zfp[l, h- J] = {Z[{k},{hi}])^^ = [ DQDQD^'^Dcp'D^pDN 



xexplzJ2jJ dt,dt2 E Qr'"'(^i,i2) [QT'"'{ti,t2) - <PT{ti)<PZ{t2. 

»2 r roo C=i,C'=i p 

- E / / dt,dt2 E n QZ ' (^1' ^2) 

"P- ° €{0,1} r=l 



il...ip 

oo 



+ / dtut)<p',{t) + thi{t)<pl{t)] 



^ fOO 

+^E<^40l(^o) -<^i] +W dtN{t) 
i=i -^0 



N 



i=l 



(43) 



We note that the segment representing stable low temperature fixed points in the (Ai, 
Aq) plane includes the point Ai = and Aq = —1. This is the only point in the segment 
that represents freezing of the correlation, a property of glasses. 



4. Time reparametrization symmetry 

We now evaluate the effect of a reparametrization t ^ s{t) of the time variable on the 
stable fixed point generating functionals. For this purpose we consider a monotonously 
increasing function with the boundary conditions s(0) = and s(oo) = oo, which 
induces the following transformations on the sources, 

Ut) = ^h{sm (44) 
kit) = hi{s{t)). (45) 
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First we consider the effective generating functional of the high temperature fixed points. 
We evaluate the fixed point generating functional of the transformed sources 

Zfp^,h]T] = j DQDQD^^D^^D~^DN 

X exp ^ E / / dt,dt2 E Qi ih, t2) (QT'^^'iti, t2) - iPriti)i^i'{t2) 



i=l ''^ 

+ / dtUt)i;^{t)+ih,{t)i;lit)] 
Jo 

+tY.^i[ipl{to)-^i]+t / dtN{t) 



N 



.1=1 



(46) 



Here wc have used new dummy variables 0, Q, Q and A^, instead of 0", 0, Q, Q 
and N, respectively, in the functional integral. We now perform the following change of 
variables 



Qr'{t,t') 



di 



' ds\" I ds 
di) [di' 



Qr'{s{t),s{t% 



Z. Oi,OC 



ds 



N{t) = ^N{s{t)), 
The change of variables results in Jacobians in the differentials, 



L)V D^)' DN 
D(t)^ Dcf)^ DN 



(47) 
(48) 

(49) 

(50) 
(51) 

(52) 

(53) 
(54) 



D0 = D(p. 

Since the field transformations are linear, the Jacobians depend only on the 
reparametrization s{t). Therefore, they are independent of the fields and sources, and 
can be taken outside the integral as common factors. 

By inserting the values of the transformed sources and dummy variables back into 
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the fixed point generating functional we obtain, 

Zfp[l,h]T\ = DQDQD(f)'^D(f)^D<fDN 



» a^,a^ \ 



dt' 



Q.-"Ks{t),s{t')) 



^0 



0/ 



^0 



+ / dt 



^iMt)msit))+Msit))^<p]{s{t)) 



N 



+i y: <Pimsm -ipi]+i dt—N{s{t)) 

^0 dt 



ds 



N 



j:ms{t))r-N 



.1=1 



(55) 



So tlien tlie transformed fixed point generating functional is 
Zfp[I,h;T] = J1J2J DQDQD(j)°D(l)^DipDN 

^^E / r^'^^'T. Qf'^'i^,^') (qT'%,s') - <Pr{s)<pfis') 



dsUs)<p',{s)+thi{s)<Pl{s)] 



N .00 
i=i -^0 



N 



.i=l 



(56) 



Here we have used the fact that a -\-a = 1. We notice that the term describing the 
coupling to the bath is not invariant with respect to the transformation t — )■ s{t), except 
in the trivial case s(t) = t. So the high temperature fixed points are not invariant under 
reparametrizations of the time variable. For the same reason, the fixed point actions 
containing both the coupling to the thermal bath and the spin-spin interaction are not 
invariant under time reparametrizations. 

Finally, we consider the fixed point generating functional for the low temperature 
fixed point family. We evaluate the fixed point generating functional for the new sources 

tJ2jJ dt,dt2 E Qi ^2) (QT'^^'iti, t2) - i^T{tiW{t2) 



-^E^n...v// «2 E nor'"M^i,^2) 



p\ 



i-i...ip 

00 



C=l,C'=l p 



e{o,i}r=i 



+ / dt[k{t)^^'',{t) + lh,{t)t^l{t)] 



N ,00 \ N ] 

+^E<^JV'- (^o) dtN{t) Y.i^l{t)f - N 

i=l -"^ li=l 



(57) 
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Here we have used the same dummy variables used in the analysis of high temperature 
fixed point actions. Perfoming the same change of variables we get the Jacobians Ji 
and J2. 

By inserting the values of the transformed sources and dummy variables back into 
the fixed point generating functional we obtain, 



Zfp[l, h; J] = J1J2 J DQDQD(l)'^D4)^D(pDN 



oii+ai 



dp I 



Qr^{s{t),s{t')) 



X QT'''^{s{t),s{t'))-ci>r{sm':^{s{t')) 



f ll...lp 

+ / dt 
Jo 



a,,,<e{0,l}r=l yO'V 



m,s{t')) 



-iMt)mm)+ihMt))-<Pi{s{t)) 



+iY,^,[<Pl{s{0)) -cp,]+i dt-N{s{t)) 



N 



Y.{<t^]{s{t))f - N 



.1=1 



(58) 



We now use the fact that a + a = 1 and that the constraints C = C = 1 ensure that 
n (if) '(^) "^"if^'t^ write the transformed fixed point generating functional 

r=l 

Zfp[l, h] J] = J1J2 J DQDQD(j)^D(f)^D<fDN 





«2 r r°o c=i,c'=i P „ , 



+ / dsUs^'^{s)+th,{s)<Pl{s)] 



roo 

+iJ2<fi[<P'M-'Pi] + i dsN{s) 
i=i Jo 

In other words, we have shown that 

Zfp[l,h; J]^JiJ2Zfp[l,h; J]. 



N 



.1=1 



(59) 



(60) 



We know that in the absence of sources, the transformation leaves the generating 
functional unchanged. This implies that J^iJ^2 — 1, but since J^i and J^2 are independent 
of the values of the sources, then for any value of the sources the fixed point generating 
functional is unchanged by the transformation, i.e.. 



Zfp[l,h] ^Zfp[l,h]. 



(61) 
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Therefore, for the Langevin dynamics of the p-spin model, the low temperature 
long-time fixed point dynamic generating functionals are symmetric under time 
reparametrizations. 

5. Discussion and conclusion 

In our long time renormalization group analysis we have shown that there are three 
families of stable fixed point dynamic generating functionals for the Langevin dynamics 
of the p-spin model: (i) a family of high temperature fixed points, which are not 
invariant under global reparametrizations of the time variable, characterized by the 
presence in the action of the coupling to the thermal bath and the absence of the spin 
interaction term; (ii) a family of low temperature fixed points with actions containing the 
spin interaction term but not the coupling to the bath, which are invariant under global 
time reparametrizations in the long time limit; and (iii) a third family of stable fixed 
points, for which both terms are present in the action, and thus the action is not invariant 
under time reparametrizations. Since not all of the stable fixed points in the model are 
invariant, it is clear that time reparametrization symmetry is a nontrivial property of 
the low temperature, interaction dominated dynamics. It should also be pointed out 
that in another interacting spin model, the 0(N) ferromagnet, the symmetry is not 
present in the asymptotic long time Langevin dynamics, even in the low temperature 
case [25]. 

The proof of invariance for the low temperature, long time dynamics only assumes 
that the couplings Ji^...ij, are uncorrelated Gaussian random variables with zero mean, 
but no condition is imposed on the variance of the couplings, thus allowing them 

to have an arbitrary space dependence. In particular, the proof applies to both short- 
range and long-range models. Since some versions of the p-spin model share many of 
the main features of structural glass phenomenology [211 [22] , we expect that analytical 
tools similar to the ones used here can uncover the presence of time reparametrization 
symmetry in models of structural glass systems. 

As discussed in Refs. [TTl UHl |19], time reparametrization symmetry is a 
spontaneously broken symmetry in a glass. The symmetry is broken by correlations 
and responses. To illustrate the spontaneous breaking of the symmetry, we consider the 
correlation function C{ti,t2)- If correlations were invariant under the transformation 
we would have C{t,t') = C {h{t) , h{t')) for all t and t' and all reparametrizations 
and the only way this is possible is when the correlation function is independent of 
time. This is not the case in glasses because the correlation decays with time. The 
presence of a broken continuous symmetry in the absence of long range interactions or 
gauge potentials is expected to give rise to Goldstone modes [27J. In the case of the 
glass problem, the Goldstone modes should be associated with smoothly varying local 
fluctuations t — )■ hr{t) in the time reparametrization UHl HH]- These fluctuations 
can be interpreted as representing local fluctuations of the age of the sample [TTl HE] ■ 
Support for this point of view comes from simulation results both in the Edwards- 
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Anderson model of spin glasses [TS] and in models of structural glasses [231 IM]- The 
kind of analysis performed in |23l [23] can in principle be straightforwardly extended to 
be applied to particle tracking experimental data showing dynamical heterogeneities in 
glassy colloidal systems [9l [10], and in granular systems [14j. 

We conclude by noting that this work hints at the possibility of analytically proving 
that time reparametrization symmetry is present in structural glasses. By investigating 
the Goldstone modes predicted as a consequence of the symmetry this may provide 
an avenue to compute detailed predictions for probability distributions and correlation 
functions that describe the behavior of dynamical heterogeneity. 
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